Supernova neutrinos, from back of the envelope to supercomputer 
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While an understanding of supernova explosions will require sophisticated large-scale simulations, it is never- 
theless possible to outline the most basic features of the neutrino emission resulting from stellar core collapse with 
a pedestrian account that, through reliance upon broadly accessible physical ideas, remains simple and largely 
self-contained. 



1. BACK OF THE ENVELOPE 

Because the central features of a system can 
often be elucidated with simple models, and be- 
cause neutrino emission is — from nature's point 
of view of raw energetics — the central feature of 
core-collapse supernovae, one may hope to un- 
derstand the basic features of supernova neutrino 
emission quickly and easilyll 

Shortly after the discovery of the neutron in the 
early 1930s Baade and Zwicky declared: "With 
all reserve we advance the view that supernovae 
represent the transitions from ordinary stars to 
neutron stars, which in their final stages consist 
of extremely closely packed neutrons" [9]. For 
the 'core-collapse' varieties of supernovae (Types 
Ib/Ic/II) this basic picture prevails today with 
strong theoretical and observational support. 

The storjH of how an 'ordinary star transi- 
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2 A self-contained but simple derivation will be carried on 
in these footnotes. In terms of nuclear physics, weak in- 
teractions, and hydrodynamics relevant to stellar collapse 
and bounce, as well as spatial neutrino diffusion during 
post-bounce delcptonization, this presentation exhibits a 
cavalier disregard for detail that borders on criminal neg- 
ligence. For responsible treatments of these aspects that 
nevertheless retain a somewhat analytic flavor (though re- 
lying to an extent on numerical simulations or experimen- 
tal information), see for instance Refs. [1 2 3 4 5 6 7fS]. 
3 In slashing through this problem, the virial theorem 10, 



tions to a neutron star,' to paraphrase Baade and 
Zwicky, involves the emission of another electri- 
cally neutral particle 'discovered' (theoretically, 
if not experimentally) in the 1930s: the neu- 
trino. Once conditions for their production are 
reached, the weakness of neutrino interactions 
implies that their emission is the most efficient 
means of cooling. Neutrino emission is responsi- 
ble for the degeneracy of the precollapse stellar 
core of mass M « 1.5 M Q (the Chandrasekhar 
mass); the core's eventual instability, at least in 
part; and the energy loss of about 1.7 x 10 53 erg 
that ultimately allows collapse to proceed to ra- 
dius -R C oid ~ 11 krn as a cold neutron star@ 

will prove to be a particularly sharp machete: 

2E N +E R = -E G , (1) 

where Ejq and En are the contributions to a star's inter- 
nal energy in nonrelativistic and ultrarelativistic particles 
respectively, and 



is the total gravitational energy of a star of uniform density 
(assumed here throughout) with mass M and radius R. 
(In this presentation G is Newton's constant; moreover, 
h = c = k = 1.) 

4 The assumption that the stellar core burns all the way 
to |gFe — the top of the binding energy curve, from which 
no further energy can be extracted by fusion — allows an 
iron core temperature Tp c 0.9 MeV to be estimated 
as follows. Fusion requires tunneling through a Coulomb 
barrier. Take the uncertainty principle Ax Ap £3 1 
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as a measure of the potency of tunneling, with Ap ft 
1jj,Z\Z^a I Ax as a momentum scale corresponding to 
the energy height of the Coulomb barrier (here /j is the 
reduced mass of colliding nuclei of proton numbers Z\ and 
Z2, and a 1/137 is the fine structure constant). Take 
also T ft Z1Z20/ 'Ax as the characteristic thermal energy 
per particle needed to bring the colliding nuclei within Ax 
of each other. Elimination of Ax yields T = 2Z' 2 Z 2 Q 2 p. 
Burning tends to proceed by successive capture of a par- 
ticles (which have modest Z, and are tightly bound and 
therefore abundant), so taking Z\ = 24, Z2 = 2, and 
p = rrtjv(4)(52)/(4 + 52) appropriate to a putative final a 
capture to ||Fe gives T Fo w 0.9 MeV (m N ft 939 MeV is 
a generic nucleon mass). 

If the electrons are not already degenerate at this point, 
they must soon become so, thanks to neutrino emission. 
A temperature Tp e ft 0.9 MeV is almost twice the elec- 
tron mass m e , so in the absence of an electron chemical 
potential of similar magnitude positrons will be thermally 
produced. But these e + will annihilate with e~ into vv 
pairs that escape freely; hence we expect cooling by vv 
emission back down towards T ft m e ft 0.5 MeV. In re- 
sponse to the resulting loss of pressure support the core 
contracts to higher density, with the loss of gravitational 
potential energy increasing the temperature and bringing 
about further pair emission, and so on. Continued increase 
in density with temperature regulated to ~ 0.5 — 1 MeV 
eventually results in electron degeneracy, which becomes 
the dominant source of pressure support. (Iron core cen- 
tral temperatures of ~ 0.5 — 1 MeV are in fact seen in 
stellar evolution codes |11|.) 

The consequences of degeneracy are dire. As the mass 
and density of the iron core increase during the final burn- 
ing stage, the velocity of the degenerate electrons ap- 
proaches the speed of light, at which point they can pro- 
vide no further pressure support. To estimate the core's 
maximum mass, take Ejy = and Br = (E F ) e in Eq. 
JTJ, where 
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is the internal energy of a uniform core supported by a 
degenerate Fermi gas of relativistic electrons; the electron 
fraction Y e is the ratio of the net number density of elec- 
trons n - — n e + (or, by charge neutrality, the number 
density of protons n p ) to the total baryon number density 
n. Note the cancellation of radius that then results in Eq. 
JTJ, a clear indication that something goes haywire when a 
star tries to support itself by relativistic degeneracy pres- 
sure alone. The resulting expression nevertheless yields 
an estimate of the mass of this extreme configuration, the 
so-called Chandrasekhar mass: 
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where Y e fi ft 0.46 is the ratio Z/A of the proton and mass 
numbers of |gFe. This may be taken as an estimate of the 
mass of the core upon collapse, and of the compact rem- 



nant that results therefrom. (Taking inhomogcncous stel- 
lar structure into account yields M ft 1.22M Q (Y e o/0.46) 2 

Next consider the role of neutrinos in the ultimate in- 
stability of the core. As it approaches the Chandrasekhar 
mass, one proximate cause of the core's instability — which 
in fact will be used here to define the onset of collapse — 
is electron capture, which saps pressure support from the 
core. Electron capture becomes possible when the electron 
Fermi energy 
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rises above the energy threshold required to convert a 
proton to a neutron. Crudely considering a nucleus of 
proton number Z and mass number A as comprising de- 
generate Fermi gases of nonrelativistic protons and neu- 
trons contained in a volume of radius r A = rzvA 1 / 3 
(where rjv ~ 1-2 fm), the typical emitted neutrino en- 
ergy l v<1 is not the typical energy of an electron captured 
from its (ambient) degenerate Fermi sea, because the en- 
ergy required to lift a bound proton from its Fermi sea 
(within the nucleus) to the top of the neutron Fermi sea 
(within the nucleus) must first be paid. Hence e l/e ft 
(e F ) e -A - [(e F ) n]A ,z ~ {£f) p; A,z] (where (e F ) n;A ,z and 
(^f)u;A,z are the Fermi energies of the neutrons and pro- 
tons bound in the nucleus), or 



(e F ) e -A- 
(e F ) e -A- 



3tt z 
I V 2 

3 / 3tt 2 



1/3 



1/3 



( A _ Z) 2/3_ Z 2/S 



m N r* N . 
(l_y e )2/3_ y<! 
m N r 2 N 



2/3 



(6) 
(7) 



where A = m n — m p — m e ft 0.8 MeV and it is assumed 
that all baryons are locked up in nuclei characterized by 
the same Z and A. Electron capture can proceed when 
the right-hand side is positive; using (e F ) e as given in Eq. 
J5J1 yields Ro ft 810 km at the onset of electron capture 
for M = 1.5 Mq and Y e ,o = 0.46, corresponding to a 
density of po ft 1.3 X 10 9 g/cm 3 . Incidentally, we are 
now in a position to compute the entropy per bayon of 
the precollapse core: sq ft 0.85 — 1.5 (corresponding to 
To ~ 0.5—1 MeV) receives contributions from the nuclei of 
mass number A = 56, the relativistic degenerate electrons, 
and photons: 
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where m A ft! m^r A is the mass of a nucleus of mass number 
A and (e F ) e = (3n 2 p Y e fi/m N y/ 3 ft 4.4 MeV is the 
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initial electron Fermi energy. (As was the case with To, 
these values of po and so are in reasonable agreement with 
the central iron core values produced by stellar evolution 
codes 

Having characterized the precollapse core and deter- 
mined that it is doomed to collapse — and supposing that 
the degeneracy pressure of nonrelativistic nucleons will 
halt collapse — what might we expect the final radius -R co ld 
of the compact remnant to be, and why do we expect it 
to consist of neutrons rather than a mixture of neutrons, 
protons, and electrons? Aware that weak interactions can 
interchange neutrons and protons, we provisionally take 
our star of "extremely closely packed [nucleons]" to con- 
sist of nonrelativistic degenerate Fermi gases of neutrons 
and protons and a relativistic degenerate Fermi gas of elec- 
trons. Hence Ejv = {Ep) n + (Ep) p in Eq. JT), where 
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and Er = (Ep) e of Eq. {3}. A second constraint is pro- 
vided by chemical ('/3' or 'weak') equilibrium, which re- 
quires 



— fl*p ~t~ /^e , 



(13) 



where the chemical potentials in include particle masses. 
That is, Hi = rrii + (ep )i at zero temperature, with 
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and (ej?)e given by Eq. J5}. Taking M = 1.5 Mq, simul- 
taneous solution of Eqs. JT) and H13I I as specified above 
yields i? co ld = 11 km and Y e co \^ = 1.0 X 10 — 2 . Hence the 
label 'neutron star' is well deserved! The single equation 
(Ep)„ = -E G /2 yields 



11 km 



(16) 



with Y e neglected. (The resulting baryon number den- 
sity n 0.25 fm -3 is within a factor of two of the value 
0.14 fm — 3 derived from from experimental measurements 
of the size of atomic nuclei, which also feature densely 
packed nucleons. A neutron star is not unlike an atomic 
nucleus of astronomical proportions!) 

Note that neutrino emission is necessary not only to get 
collapse started, but also for it to proceed to completion. 
This is another consequence of the virial theorem: from 
Eq. (TJ we see that the total energy of the precollapse 
core is (Ep) e + Eq 0, while the total energy of the cold 
neutron star is (Ep) n + Eq Rj Eq/2. Hence an amount 
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Figure 1 . In the crude model of gravitational collapse 
presented here, the neutrino diffusion time scale t„ 
exceeds the free-fall time scale tg when the core ra- 
dius has shrunk to about 180 km, which corresponds 
to a trapped lepton fraction of Ye, trap « 0.33. 

Basic characterizations of the precollapse core 
and the cold neutron star to which it collapses 
may be sufficient to estimate the total energy 
loss to neutrinos, but additional features of this 
emission — which flavors are emitted, at what en- 
ergies, over what time scales — depend crucially 
on the flow of electron lepton number and the gen- 
eration of entropy. 'Neutronization' via electron 
capture entails the flow of electron lepton number 
from degenerate electrons to electron neutrinos, 
and the thermal energy associated with finite en- 
tropy results in the thermal emission of neutrinos 
and antineutrinos of all flavors @ 

The salient point about the flow of electron 
number during collapse is that while electron neu- 
trinos initially escape freely, they eventually be- 
come trapped, as Fig. Q] illustrates^ A roughly 



of energy equal in magnitude to half the final neutron star 
gravitational energy, 

„ / M \ 2 / 11 km\ 

E ^ 1 - 7xw ^(nwJ (— )' (17) 

must be lost to neutrinos (and a bit of it to the explosion) 
before the neutron star can reach its cold final state. 

5 In a sense these two themes are complementary, for de- 
generacy entails low entropy in general, and in particular 
a net electron number prevents the rise of a positron pop- 
ulation that would allow vv pair emission of all flavors by 
e+e — annihilation. 

6 That neutrino trapping will occur is shown by a com- 
parison of the collapse and neutrino diffusion time scales, 
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Figure 2. Crudely estimated neutrino luminosities 
(thick) and characteristic energies (thin). Note the 
change in time and luminosity (but not energy) scales 
between the left and right panels. The left panel is a 
close-up of infall and bounce at t = 0. 



but before examining these a word on weak interactions is 
in order. For energies much lower than the masses of the 
W + '~ and Z° bosons (whose exchange are responsible for 
'charged-' and 'neutral-current' interactions respectively), 
the effective weak interaction Lagrangian is 

C = -(4G F /V / 2)(Jcc)j J (Jcc) M - (G'F/v / 2)(Jnc) M (Jnc) M , 

where the Fermi constant Gp = 1.17 X 10 -11 MeV -2 sets 
the overall scale of weak interaction amplitudes (/3 decay 
being the prototype), and (J CC ) M and (Jnc) 11 are a 'charged 
current' (which changes electric charge by one unit) and 
'neutral current' constructed of quark and lepton fields 
(e.g. Appendix B of Ref. |13|). In this presentation the 
factors SG 2 F and G F /2 are taken as crude estimates of 
leading factors in charged and neutral current processes 
respectively, with neutrino energy factors supplied to give 
the correct dimensionality for rates or cross sections. 

Now to the collapse and neutrino diffusion time scales. 
Assuming that the collapse attending the core's approach 
to the Chandrasekhar mass is catastrophic, such that the 
core is essentially in free-fall, the instantaneous contrac- 
tion time scale when the core is at density p is 
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On the other hand the trapped electron neutrino diffusion 
time scale is 
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where \ Ue is the neutrino mean free path, o v a—>Av ~ 
A 2 (G F /2) e 2 ,^ is the cross section for coherent neutral 
current scattering off a heavy nucleus of mass number A 
(the dominant source of opacity; the argument that the 



estimated 5.0 x 10 50 erg are released in v e before 
trapping; the gradual rise in total luminosity and 
cutoff at trapping, as well as the typical neutrino 
energy during infall, are depicted in the left panel 
of Fig. [2] at negative (i.e. pre-bounce) times0 



baryons remain locked in heavy nuclei through trapping 
is given in the next footnote). In accordance with the 
above remarks, a neutral current scattering cross section 
per nucleon (G%/2) ej^ is suggested by dimensional anal- 
ysis. Applied to an entire nucleus yields not just a factor 
A corresponding to an incoherent sum over scatterings off 
individual nucleons, but a factor A 2 corresponding to col- 
lective scattering off the entire collection of nucleons — a 
coherent sum in the scattering amplitude. This is be- 
cause the neutrino wavelength is larger than the size of 
the nucleus (as can easily be verified from Eq. (0 and 
TA = rjvA 1 / 3 with rjv ~ 1.2 fm). With the nucleus then 
treated as a point particle, the cross section is partially 
analogous to Rutherford scattering of a charged particle 
from a point nucleus of proton number Z, with our factor 
of A 2 associated with 'weak charge number' A correspond- 
ing to the factor Z 2 in the Rutherford scattering formula. 
As collapse begins, 
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and r Vc K, since e„ e as given by Eq. Q has here 
been taken to define the onset of collapse. While neutrinos 
initially escape easily, the decrease in tg and increase in 
t Vc as R decreases during collapse imply that at some 
point neutrinos will be trapped (in addition to the explicit 
dependence of T„ e on R in Eq. il9l , see Eq. (0 for the 
increase in (ep) e and hence e v<! as R decreases). 

In order to estimate the trapped lepton fraction, note 
that as long as the u e from electron capture escape freely, 
the electron fraction changes according to 
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The factor of Y e gives the fraction of baryons that are 
protons; the charged current capture rate per proton on 
nuclei of mass number A, r eA —>Av ~ 8G|,e®, has been 
estimated in accordance with previous remarks (in con- 
strast to neutrino scattering, there is no coherence on the 
nucleus as a whole) ; and dt — drQ has been taken. Equa- 
tion H22I I. with Eqs. (0 and J5}, is an ordinary differential 
equation for Y e whose solution is shown in Fig. [l] Using 
this solution in connection with Eqs. J 18D and 11191 . the 
trapped lepton fraction is the value of Y e where tg = t v . 
This is illustrated graphically in Fig. [T]for M = 1.5 Mq, 
R = 810 km, and A = 56. The result is y e ,tra P « 0.33. 
7 The total energy released in u e before trapping is ob- 
tained by integrating dE Ve = e Vs (M /mi^)(dY e /dR)dR 
from the initial radius to the trapping radius. The lu- 
minosity during infall is L„ c = dE V(1 /dTG, and is plotted 
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The other point to make about collapse is 
that the entropy per baryon changes only mod- 
estly, which conditions what happens on dynam- 
ical time scales at the end of implosion|f| Instead 

in Fig. |2|as a function of time t = —tq. The neutrino en- 
ergy is given by Eq. J7jl . In all of this the solution Y e (R) 
as obtained in the previous footnote is used. 
8 To see that the entropy per baryon changes only modestly 
during collapse, consider the first law of thermodynamics 
in the form 

Tds = dq - fiidYi. (23) 

i 

Here dq = d(e /n) +p d(l / n) (where e is the internal energy 
density, n is the baryon number density, and p is the pres- 
sure); but it is equal, as the label dq suggests, to the heat 
per baryon exchanged with the environment. The second 
term represents the entropy per baryon generated by out- 
of-equilibrium 'chemical' changes to the core (fi{ are the 
chemical potentials including rest masses of the species 
i, and dYi are the changes in the numbers per baryon). 
There are arguments for the smallness of ds both before 
and after trapping. 

Before neutrino trapping both energy and lepton num- 
ber are lost via electron capture, and the modesty of 
entropy increase results from a substantial cancellation 
between these contributions. The energy per baryon 
lost to electron capture neutrinos that freely escape is 
dq ss t Uc dY e . Meanwhile the chemical change of elec- 
tron capture involves dY e = dY p = — dY n (the first 
equality arising from charge conservation and the second 
from baryon number conservation), so that y\ [ijdYj = 
(p e + (i p — [i n )dY e . Hence 

Tds rs (e„ e + [i n - Hp — [i e )dY e . (24) 

It turns out that the leading factor cancels exactly with 
the prescription for electron capture neutrino energy of 
Eq. 10. In a more careful calculation the cancellation 
would not be exact, but it would still be substantial. An 
important consequence of the small change in entropy be- 
fore trapping is that the baryons remain in nuclei (which 
feature the coherent scattering cross section oc A 2 ): using 
A = 1, R = 180 km, and T = 1 - 5 MeV in Eq. © 
shows that s R3 4 — 6 would be required for the nuclei to 
disassemble into free nucleons before trapping. 

It can also be argued that the change in s between trap- 
ping and core bounce will be modest, this time because 
the terms of Eq. I|23| l tend to vanish individually. First, 
dq ss since neutrinos can no longer carry away energy. 
Second, with neutrino trapping, 

[li.dYi = (/i e +fip- [i n - [iv^)dY R + [i Vc dY L , (25) 

i 

where Yi = Y e + Y„ E is the trapped lepton fraction (in 
analogy with the definition of Y e given earlier, Y Ve = 
(n Ve — na e )/n). Between trapping and bounce, while 



of a soft landing cushioned by abundant degrees 
of freedom capable of gently absorbing the ki- 
netic energy of infall, a hard 'bounce' can be ex- 
pected when nucleon degeneracy pressure kicks in 
at supranuclear densities]^ This is accompanied 
by a burst of v e emission (see the left panel of Fig. 
[2]), involving an energy loss of about 5.9 x 10 51 erg 
and the sudden release of trapped electron lepton 
number F°l Moreover, an additional 6.8 x 10 51 erg 

Yi, = y ej trap, the second term vanishes because dYi Ri 0, 
and the first term also becomes small as e~ , e+, p, n, v R 
and P e come into chemical ('/?' or 'weak') equilibrium. 

9 The estimates of post-bounce neutrino emission pre- 
sented below involve reasoning from the approximate 
global states of a neutron star in chemical equilibrium at 
finite temperature. In addition to thermal photons and 
Vy.Vn and v T v T pairs, there must be allowance for arbi- 
trary degeneracy in y e P e pairs, relativistic e + e~ pairs, 
and nonrelativistic neutrons and protons. For relativis- 
tic pairs of arbitrary degeneracy there are simple ex- 
pressions for the differences of number densities and the 
sums of energy densities of particles and antiparticles, as 
well as the entropy 6 3. For the nonrelativistic nucle- 
ons the needed quantities can be derived from the Lan- 
dau potential Q(T,V, /J,) = -pV = E -TS - [iN . This 
can be written in terms of a polylogarithm function as 
fi(T, V,[i) = V(m 3 T 5 /2 7 r 3 ) 1 /2Li 5/2 [- expf>/T)], whose 
symbolic derivatives and numerical evaluation can be han- 
dled by any self-respecting symbolic mathematics package 
(though the Sommerfeld expansion seems more reliable for 
[i/T > 10 — 20). The global neutron star states are charac- 
terized by two parameters (the entropy per baryon s and 
the lepton fraction Yi = Y e + Y Ve ) and seven unknowns 
(R, T, y e , fi e , [i Ve , [J-n, and [i v ). The seven equations to 
be solved are the virial theorem of Eq. JT} ; chemical equi- 
librium, Eq. H13I I with [i Ue added to the left-hand side; the 
total entropy per baryon; and four 'particle number equa- 
tions' defining Y e , Y Vc = Yi,— Y e , Y p = Y e , and Y n = l—Y e 
in terms of the chemical potentials, temperature, etc. 

10 A neutron star state (as described in the previous foot- 
note) with s = 1.5 and Yi = Y e ,trap has a total energy 
Etot much lower than -£„ e) infall (the energy lost to u e es- 
cape during infall), so it cannot be the post-bounce state. 
(As noted earlier, the precollapse initial condition sup- 
ported by purely relativistic energy is characterized by 
Etot = 0, which follows from Eq. lHJ.) Examining a 
sequence of states with higher s in a bid to reach Etot = 
— E v i n f a n eventually yields a 'burst state' with s ss 5.9 
whose radius R a; 303 km is equal to its electron neutrino 
trapping radius. (As before, the neutrino trapping radii 
are given by equating a diffusion time scale r u = 3R 2 /X 
with a dynamical time scale Td yn = (Gp) -1 / 2 . Here 
A„ e ~ (njvcrjv +n„<j n )~ 1 , where crjv ~ Gf€ 2 c /2 is for 
neutral current scattering on nucleons and <j n Ri 8Gfe 2 
is for charged current absorption on neutrons. For later 
use, note that Xp e as (njv""JV +npff P ) _1 , where <r p cr n , 
and X u & u T ,o r ~ (™jv°"jv) _1 -) Hence as the core re- 
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or so is lost from the core before the dust set- 
tles after a dynamical time scale of about 24 ms 
post-bounce! 11 ! in the left panel of Fig. [5] this full 
amount is attributed to neutrino emission equally 
divided among all flavors f 1 ^! 

The neutrino emission encountered thus far ac- 
counts for only about 10% of the total energy to 
ultimately be lost, with the rest — about 1.5 x 10 53 
erg — temporarily stored as heat, thanks to the en- 
tropy generated by bounce and its aftermath. As 
the hot nascent neutron star shrinks (see Fig. [3]), 
this energy is lost to neutrino emission of all fla- 
vors (albeit with a modest hierarchy of luminosi- 
ties and average energies) over a period of many 
seconds, as depicted in Fig. [2] — a recognizable if 
very rough caricature of light curves produced by 
detailed models I14I15I16I17I P^l 



bounds to this radius, u e with Si 3/i e /4 ~ 7.1 MeV 
escape freely in a burst in which a lepton fraction dY^ 
and energy E u burst ~ £i> e {M/mtf)dYi, are lost during 
a time 7d yn Si 24 ms that characterizes this 'burst state.' 
In order to find dY^, a sequence of states is searched by 
decreasing Yl while adjusting s to stay just below the 
(varying) v e trapping radius. All such states are found to 
be characterized by Etot < — -Ev c ,i n fall — E„ e burst > aii the 
way down to Y Ve ~ 0. The conclusion is that in this crude 
model the deleptonization accompanying the v e burst is 
essentially complete, with E^burst ~ (-^/ m iv)^e,trap 
evaluated at the initial 'burst state.' 

11 At the end of the sequence of states with declining lep- 
ton number searched as described in the preceding foot- 
note, we still do not have a state with -Etot = — E„ ei i nla n — 
E„ e burst; further energy loss is needed. But it is found 
at the end of this searched sequence that the population 
of thermal neutrinos of all flavors has begun to make a 
nontrivial contribution to Etot, so their emission is now 
a viable means of prompt energy loss so long as the star 
is outside the trapping radius of these flavors. Hence our 
'post-bounce state' after the dust settles is that character- 
ized by Yjj = Y e and the highest entropy per baryon that 
keeps the star within the smallest trapping radius (that 
of Vfj. t T and their antiparticles, which as described in the 
previous footnote have the largest mean free path). This 
is the state with s si 5.3 and R si 130 km at the far right 
of Fig. |3] The additional energy loss referred to in the 
main text is the difference between Etot of this state and 

E[y ei i n f a n Ejy e burst • 

12 Some fraction of this would actually be transferred on 
a dynamical time scale to layers exterior to the core — in 
an incipient 'prompt' explosion, it was (rashly) hoped in 
bygone days — by processes that cannot be discussed in 
detail here, such as the direct momentum impulse of the 
rebounding core, 'pdV work, shock heating, and so on. 

13 The sequence of neutron star states in Fig. [3] is that 
characterized by Yl = Y e and s stepped down from 
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Figure 3. As the nascent neutron star shrinks in 
radius in response to neutrino emission during the 
post-bounce cooling phase, the total energy _Etot of 
the core and entropy per baryon s decrease mono- 
tonically even as the temperature T increases before 
making its final plunge. Starting at 0.1 s post-bounce, 
diamonds mark time intervals of 0.2 s; starting at 1 
s, squares mark time intervals of 2 s. 

2. SUPERCOMPUTER 

That some basic features of supernova neu- 
trino emission can be motivated by a stick-figure- 
quality model provides a baseline degree of in- 
stant gratification. This is in contrast to the ex- 
plosion manifest in an expanding supernova rem- 
nant's kinetic energy, which, as a 1% subsidiary 
detail from nature's perspective, has been requir- 
ing decades of detailed study by numerous people 
to understand 18 



5.3 to 0. The luminosities L Vi and average energies 
for the various species are found by equating two differ- 
ent expressions for the neutrino luminosities necessary to 
allow a transition between two successive neutron star 
states. Defining dEtot to be the difference in Etot be- 
tween two successive states in the sequence, and assum- 
ing that the star contrives to release an equal amount 
of energy in each flavor, one expression is the ratio of 
— dEtot/6 to the neutrino diffusion time scale r v (given for 
the various species in a previous footnote) at that point 
in the sequence. A second expression for luminosity is 
L„ 4 = 47ri? 2 (77r 2 /960)(?„ i /3) 4 , a neutrino version of the 
Stefan-Boltzmann law with T e ff j ss e lJi /3. Once the e Vi 
and L Vi are determined, the sequence is related to a time 
coordinate by defining dt = — dEtot / (^ . L Uj ) for each 
step of the sequence. 

14 Such is the (probably unavoidable) result of our anthro- 
pocentric chauvinism. For decades, and not yet sure even 



7 



But stick figures do not ultimately satisfy, of 
course; we long for the voluptuous contours of 
a fully- fleshed-out model. In fact, under thor- 
ough dissection the simple model of the previous 
section might stand revealed as more retrospec- 
tively illustrative than predictively robust 13 This 
is one indication of the ultimate necessity of de- 
tailed simulations. Another is that the details 
of early neutrino emission and the explosion are 
more interrelated than some of my previous com- 
ments might suggest. Whether it be by playing 
the spoiler, acting as agent of explosion, or setting 
the conditions necessary to the operation of some 
other mechanism, neutrinos have played direct or 
indirect roles in the many explosion scenarios dis- 
cussed over the years — and detailed simulation of 
collapse and the second or so after bounce will 
allow a neutrino signal detected from a Galac- 
tic supernova to serve as an unrivaled diagnostic 
tool. In its full glory neutrino transport is a time- 
dependent six-dimensional problem, as it requires 
the tracking of neutrino energy and angle distri- 
butions at every point in space. The computa- 
tional resources necessary to solve the daunting 
fulness of this problem are still just over the hori- 
zon, but efforts to meet this problem in its full 
measure are already underway [19120] . 
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